APS/ 

Characterizing heralded single photons with an all-fiber source of 

photon pairs 

Lei Yang, Xiaoxin Ma, Xueshi Guo, Liang Cui and Xiaoying l|*] 
College of Precision Instrument and Opto- electronics Engineering, Tianjin University, 
Key Laboratory of Optoelectronics Information Science and Technology, 
Ministry of Education, Tianjin, 300072, P. R. China 

Abstract 

Based on the signal and idler photon pairs produced in a piece of high nonlinear fiber by a pulsed 
pump, we characterize the heralded single photon source from both the theoretical and experimen- 
tal aspects. In the theory model, started from the derived expression of Bogoliubov transformation 

(2) 

for a broadband pulsed pump, the second-order intensity correlation function g c (0), heralding 
efficiency H, indistinguishability and brightness of the heralded single photons as a function of 
source parameters are analyzed and discussed. In the experiments, using several kinds of combi- 
nations of the source parameters, the values of gc (0) and H are measured and compared. The 
experimental results are consistent with the theoretical predictions. The investigations are use- 
ful for optimizing the parameters and for developing a single photon source suitable for quantum 
information processing. 
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I. INTRODUCTION 



; that a single isolated quantum system 
4^6\. The second uses the spontaneous 



Generation of single photons is not only important for the fundamental test of quantum 
mechanics, but also crucial for quantum metrology and quantum information processing 
An ideal single photon source would produce completely characterized single photons on 
demand. However, all the available sources fall short of this ideal. To evaluate the quality 
of a single photon source, the following three criteria are often taken into account: (i) pho- 
ton statistics, represented by the second-order intensity correlation function; (ii) efficiency, 
characterized by the probability of delivering a photon in response to a request; and (iii) 
indistinguishability, described by its mode structure. In practice, the requirements depend 
on its specific applications. For example, for linear optics quantum computing involving in- 
terference between photons from multiple sources [2j, the indistinguishability is important; 
while for quantum key distribution (QKD), which works well with a few modes, the main 
concerns are the photon statistics and efficiency 

In general, single photon sources can be classified into two categories: deterministic source 
and heralded source. The first type relies on the fac ; 
can emit only one photon at each time it is excited 
parametric emission of distinguishable signal and idler photon pairs, the detection of the 
one member of the pair indicates the existence of its twin photon. Since the physical process 
of parametric emission is constrained by the energy and momentum conservations, the de- 
tected location of the heralding photon defines the location of its twin. This is a significant 
advantage over the other type of source. 

The heralded single photon sources (HSPSs), which are inherently probabilistic, have 
traditionally been based on spontaneous parametric down conversion (SPDC) in x crys- 
tals. Since the first conditional preparation of single photon state reported by Hong and 
Mandel this kind of source has been characterized from various aspects, and signifi- 
cant progress has been achieved in improving collection efficiency, in optimizing the photon 
statistics, and in engineering the mode structure [8|. Recently, experiments and theoretical 
investigations have shown that photon pairs generated in optical fiber by y^-based spon- 
taneous four wave mixing (SFWM) can also be used to develop the HSPS|9Nl2|. Moreover, 
the fiber based HSPS has the advantages in single spatial mode, low collection loss, compact 
physical size, low pump power and compatibility with an optical network. 
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The spontaneous parametric emissions can be realized by using either a continuous wave 
laser or a mode-locked pulsed laser. For the quantum information processing applications 
requiring synchronization , photon pais with creation time defined by the pulsed 

pump are often preferred. The HSPS based on the pulsed SPDC has been experimentally 
demonstrated for some time, but its comprehensive theoretical description is not available 



until recently [14Nl7l|. In Ref. 



161 ]. the brightness and purity of HSPS are analyze d by 



taking both the spectral and spatial degree of freedoms into account; while in Ref. |17| . 
the dependence of the three criteria for evaluating the quality of HSPS are simulated by 
expanding the output state of the SPDC in x crystal based waveguide to second order 
in photon number. In contrast, for the the HSPS based on the SFWM in fibers, which is 
usually pumped with a pulsed laser, there has been a few experimental demonstrations 0- 



11] , however, the detailed characterizations have not been done yet. 

In this paper, using the signal and idler photon pairs generated in 20-meter long high 
nonlinear fiber (HNLF) by a pulsed pump in 1550 nm band via SFWM, we carefully char- 
acterize the heralded single photon source from the theoretical and experimental aspects. 
In the theory model, to take the multi-photon emission into account in a simpler way, the 
evolution of signal and idler fields is calculated in Heisenberg picture. Based on the derived 
expression of the Bogoliubov transformation for a broadband pulsed pump, the fomular 
of the second-order intensity correlation function g c (0) and heralding efficiency H of the 
HSPS are obtained with analytic method for the first time. In the experiments, exploiting 
the advantaged conditions provided by the various kinds of commercially available off-the- 

(2) 

shelf fiber components, the values of g c (0) and H are measured and analyzed by using an 
all fiber HSPS, and the results agree with theoretical predictions. 

The rest of the paper is organized as follows. Section 2 gives the conceptual representation 
of the HNLF-based HSPS. Section 3 works out a theoretical frame suitable for describing the 
HSPS, the dependence of the photon statistics, heralding efficiency, indistinguishability and 
brightness upon the source parameters are theoretically analyzed and discussed. Section 
4 demonstrates the experiments of the all fiber HSPS with several kinds of combination 
of source parameters, and the experimental results are compared with the theoretically 
calculated results with no fitting parameters. Finally, we give a brief conclusion and discuss 
the possible improvements. 
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II. CONCEPTUAL REPRESENTATION OF THE HERALDED SINGLE PHO- 
TON SOURCE 

The conceptual representation of the scheme for generating heralded single photons is 
shown in Fig. 1. The HSPS is based on the photon pairs generated in HNLF via SFWM 
by a pulsed pump with a pulse width and repetition rate of about 10 picoseconds and 41 
MHz, respectively. When the central wavelength of the pump pulses is in the anomalous- 
dispersion regime, the parametric process is phase-matched and the probability of SFWM 
is significantly enhanced. In this process, two pump photons at frequency u p \ and u p2 , 
respectively, scatter through the Kerr (x^) nonlinearity of the fiber to create energy-time 
entangled signal and idler photons at frequencies u s and u>i, respectively, such that u p i+u p2 = 
u s + 10 i. The photons created in pairs are originated from the vacuum fields at frequencies 
u s and Ui, ai n (u s ) and ai n (ui), respectively. At the output port of the HNLF, photon pairs 
predominantly co-polarized with the pump photons are described by field operators a out (u s ) 
and a out {ui). 

To measure the signal and idler photons, one must effectively suppress the pump photons 
from reaching the detector, so the output of HNLF is passed through a dual band filter F2. 
Signal and idler photons, described by the field operators, b{u s ) and b(ui), are then detected 
by the photon counting system consisting of three InGaAs/InP avalanche photodiode-based 
single photon detectors (SPDs). The SPDs are operated in the gated Geiger-mode, the gate 
pulses arrive at a rate of about 2.58 MHz, which is 1/16 of the repetition rate of the pump 
pulses, and the dead time of the gate is set to be 10 /is. The pulse widths of the gate for 
the SPDs are 2.5 ns. The timing of the gate pulses are adjusted by a digital delay generator 
to coincide with the arrival of signal and idler photons. 

The photons in idler band are detected by SPD1, whose detection signals are used to 
herald the existence of single photons in signal band. To characterize the HSPS, photons 
in signal band propagate through a 50/50 fiber coupler and then detected by SPD2 and 
SPD3, respectively. To measure the second order coherence function and heralding efficiency, 
besides recording the single counts of each SPD, the two-fold and three-fold coincidence 
between SPD1, SPD2, and SPD3 also need to be recorded. 

It is worth noting that in our experiment, the detuning between signal (idler) and pump 
photons is about 0.8 THz. In this case, the photons in signal and idler bands originated 
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from self phase modulation are negligible. However, the SFWM in HNLF is inevitably 
accompanied by Raman scattering [18]. For simplicity, only the SFWM process is considered 
in our theoretical analysis. 

III. THEORY 

To characterize the most important feature of a HSPS - - the second order intensity 
correlation function, the theoretical model should take the multi-photon emission of the 



source into account. This issue has been investigated in Ref.[19|, in which the HSPS relies on 
continuous-wave SPDC in x^ crystal, and the analysis of second-order intensity correlation 
properties is based on the signal and idler fields evolved via the Bogoliubov transformation 



in single mode form 20] . To work out a theoretical frame suitable for describing our fiber 
based HSPS, an expression of the Bogoliubov transformation for broadband pulsed pump is 
derived, and the theory in Ref.[19[ is expanded to a multi-mode form. 

;he volume within the fiber that is capable of describing 



The interaction Hamiltonian in 
the interaction of SFWM is 



21 



22] 



Hi = a X (3) J dV (E+E+E^ + h.c.) , (1) 

where a is a constant determined by experimental details, E^ (i = 1, 2) is the positive 
frequency electric-field operator of pumps, E~ and E~ are the negative frequency electric- 
field operators of signal and idler fields, respectively, and h.c. stands for Hermitian conjugate. 
The integral is taken over the entire volume of interaction V. Since we assume all the optical 
fields are linearly co-polarized, the tensorial nature of the nonlinear electric susceptibility 
can be ignored. 

In the model, the strong pump pulse propagating along the fiber (denoted as z direction) 
with a Gaussian shaped spectrum remains classical, it can be written as 

E+ = E p0 e-^ PpZ J du#e-^-'* o)a/2a $e ik * ig - iu * it , (2) 

where P p = 2y/nA e ffE cn(u)<T p l E p l is the peak power of the pump pulse, which is treated 
as a constant under the undepleted pump approximation, A e ff and e denote the effective 
mode area of the optical fiber and the vacuum permittivity, respectively; c and n(u) are 
the speed of light in vacuum and refractive index of the fiber, respectively; a p and co p o are 
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the optical bandwidth and the central frequency of the pump pulses, k pi is the wave vector 
of pump, and 7 = 3u p oRe(x^) / ('in( y uj) 2 c 2 eoA e ff) is the nonlinear coefficient. The signal 
and idler fields are quantized electromagnetic fields as given by the following multi-mode 
expansion: 



where a*(cjj) is the creation operator of the field at frequency Uj, kj is its wave vector, and 
Vq defines the quantization volume. 

Substituting Eqs. (EJ)- ((3J) into Eq. (TT]), and carrying out the integration over the volume, 
the interaction Hamiltonian can be written as: 

Hi = — - p [ du p idu P 2du s du i a Ji (us)^ (ursine ( J e 2 

°p J V 2 / 

X g-{( tJ pi- a 'po) 2 +( w p2-Wpo) 2 }/2o-| e -i(a;pi+tjp2-w s -Wi)t _|_ ^ c ^ 

where a' = aA e ffhc/ (3y/7reonVQ), L is the length of the optical fiber, and Ak = k pl + 
k P 2 — k s — ki — 2 r )P p is the phase-mismatching term. For the pulsed pump with the central 
wavelength in the anomalous-dispersion regime of HNLF, the bandwidth of the sine function 
is usually greater than 10 nm, which is 10 times wider than that of filter F2 used in the 
experiments. Therefore, we can assume sine (^nO ~ 1. 

To obtain the unitary evolution operator U = exp{^ J dtH^, we need to carry out 
the integration over time and over all the possible combinations of u p i and u P 2 within the 
pump bandwidth. Since the integration over time gives rise to the 5 function in Au = 
u p i + ujp2 — 0J s — (jJi to guarantee the energy conservation at single photon level, the evolution 
operator can be expressed as: 

U = exp | — J doj s duji(j)(uj s , Ui)a}{iO s )a}{uJi) + h.c. | (5) 

with 

<f>{u) a ,u>i) = exp{--^(w s + u { - 2u p0 ) 2 }, (6) 

4<T p 

where G = —i2ir<\/ : Ka!jP p L/h, is proportional to the gain of SFWM. 

In the Heisenberg picture, the field operator of signal (idler) beam at the output of HNLF 
can be written as 



L out\ 



(u s{i) ) = U ] a in (uj s{i) )U = / du' s( i ) h 1 (u' s{£) ,uj sii) )a in (uj' s{i) )+ / du i{s) h 2 (u i{s) ,u s(l) )a\ n (u i(s) l7) 



with 



1, 



(8) 



where the complex functions h^cj',^, cj s (i)) and h,2(cUi( s ), u s (i)) are given by 



»(<) 
1 



|G| 



2/i 



exp{- 



W s(i)) 2 } 



(9) 



and 



|G| 2n G 



^ V2n + l(2n + 1)! 2V7TO-p 



exp{- 



4cj2 x (2ra + 1) 



The field operators ow(u; s ) and a out {uji) commute with each other, and individually obey 
the commutation relation [a out (uj), al ut ((jj')] = 8 (to — u/). Therefore, the transformation rep- 
resented by Eqs.(17|l- (110p is equivalent to the Bogoliubov transformation. To our knowledge 
this is the first formulation of such a general expression for broadband pulsed pump. 

In the low gain regime, which is of interest to us, a out (uj s ^) can be expressed in the 
following form: 

G f 1 

aout(u s {i)) = flm(w s (i)) H / dui( s ) exp{- — (u i{s) + u s{i) - 2u p0 ) 2 }a\ n (uji {s) ) + o(G).(ll) 

a p J 4cr p 

After passing through the dual band filter F2, which is described by the function fj(oOj) = 
} (j=s,i), where Ujo and Oj respectively denote its central frequency and 



(Ujj-Uljo)' 



exp{- 

bandwidth, the field operators of signal and idler photons involve into 

b(uj) = y/rfifjiw^aiuj), (12) 

where rjj (j=s,i) is transmission efficiency of F2. 

Using the field operators described by Eqs (TTTT) and (TT2"|) . it is straightforward to obtain 
the electric-field operators at the SPDs, with which the photon statistics and heralding 
efficiency of the HSPS can be calculated. If the fcth pump pulse is centered at tofc in time 
domain, the idler field detected by SPD1 (per pulse) can be written as 



E 1 (t — t k) 



^7T 



dujib(uj.;)e 



-iu]i(t—t ()k ) 



(13) 



where 771 is quantum efficiency of SPD1; and the signal fields, passed through the 50/50 
beam splitter and detected by SPD2 (SPD3), can be written as 



E 2(3) (* ~~ 



OA.- J 



du s b(u s )e- luJs{t ~ tok \ 



(14) 



where 772(3) is quantum efficiency of SPD2 (SPD3). 

Because the individual signal and idler fields generated by the SFWM parametric process 
are thermal states, the joint state of signal and idler fields, originated from the kth and hh 
pump pulse, is a zero-mean Gaussian state. Therefore, the second order moments of the 
joint state are given by its temporal auto- and cross-correlation functions as follows: 



(0\Er( tl -t 0k )Ef(t 2 -t 0l ) |0> ex \G\ 2 J 

x / dwi( 8 )(l>*(ui(B),U8(i))<f>(ui(s),Us(i)) U = 1,2,3), (15) 



(0| Et(ti - t 0k )E+ 3) (t 2 - t ol ) \0}<xGj ^^/.(^/iCwOe^'^-^^^-*")^^,^^) 
and 

(0| E^{t x - t 0k )E+ 3) {h - hi) |0) = (0| E+(ti - t 0k )Ef{t 2 - toi) |0) = 0. (17) 
A. Photon statistics 

In the process of analyzing the second-order intensity correlation of the heralded single 
photons, we first find out the count probability of SPD1, SPD2, and SPD3, respectively; 
then, to determine the reliability of HSPS, we calculate the two-fold coincidence between 
SPD1 and SPD2 (or SPD3) to quantify the quantum correlation of photon pairs; finally, to 
show the capability of creating one and only one photon per heralding events, we compute 
intensity correlation function for the signal field, conditioned on observing an idler photon 
count. 

According to Eqs (H3&-CHI), the count probability per pulse of SPD1 and SPD2 (SPD3) 
can be written as: 

rT r /2 
r/2 



T /2 

Pi= dt (0| Ei(t - t ok )E+(t - t ofc ) |0> = V2n \G\ 2 (18) 

J-T r /2 a v 



and 

cT r /2 



/I r /Z 
dt <0| E-{t - t 0k )E+(t - t ok ) |0) = ^= \G\ 2 r)sV2(3)-, (19) 
-T r /2 V2 ff p 

where T r is the response time of SPD, determined by its corresponding gate width. 
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We note that the count probability of any SPDs is a time integral from — T r /2 to T r /2. 
Because the gate width of SPDs is much longer than the creation time period of photon 
pairs confined within the pulse duration, Pi and P 2 (3) i n Eqs. (ITHjl and (IT^|) can be treated 
as a integral from — oo to oo. 

The two-fold coincidence count probability per pulse between SPD1 and SPD2 (SPD3) 
can be expressed as 

p+Tr/2 r+T r /2 

Pi2(3)(tok,toi) = / dt s dti (0| EZ 3 \(t s — t ok )Ei(ti — toi)E 1 (ti — t i)Et 3) (t s — t ok ) |0$20) 

J -T r /2 J-T r /2 

We carry out the calculation by using the quantum form of the Gaussian moment-factoring 
theorem. For the signal and idler photons generated from the different time slots of pump 
pulses, we have t ok — t i ^ and t ok —t i ^> T r . In this case, the coincidence count probability 
is 

^12(3) (Wo*) = Pi P2, (21) 

which is originated from the multi-photon events, and is referred as the accidental coinci- 
dence. While for the signal and idler photons produced in the same time slot, i.e., t ok = t i, 
we arrive at 

Pl2(3)(0) = PlP 2 + lvsV2(3)PlL, (22) 

where £ s = , °" s = is the collection efficiency of the photon pairs for the heralding 
photons detected in idler channel [23|, and in the bracket stands for i fc —toi = 0. Eq. (122"]) 
indicates the coincidence of signal and idler photons produced by the same pump pulse is 
the sum of the accidental coincidence and the true coincidences come from the pair events. 

In carrying out the calculation in Eq. f[2"2"j) . the result of the time integral in Eq. ( 12"U|) 
becomes a 5-function. Moreover, considering the quantum correlation only exist between 
the signal and idler photons produced in the same pulse, for the correlation hereinafter 
calculated and verified, we let t k = toi = 0, and treat the time integral from — T r /2 to T r /2 
as that from — oo to oo. 

With the expressions of the two-fold coincidence count probability, we can analyze the 

(2) 

quantum correlation of the photon pairs, g si (ti , t 2 ) , which is determined by the degree of 
second-order intensity correlation between the signal and the idler fields 

(2), v ( E 2(3) ( t 2)Ei(t 1 )E+(t 1 )E+ 3) (t2)) 

9si [1,2) ~ (^r(ti)^ + (ti))(^-(3)fe)^ + (3) (t2))- 1 j 



Eq. f l23|) shows that the precise knowledge of g% (ti,t 2 ) needs SPDs with a resolution time 
shorter than the duration of pump pulse. However, in practice, limited by the resolution time 
T r , the measurement of g^(ti, t 2 ) is a time integral from —T r /2 to T r /2, i.e. (g^\ti, t 2 ) 



p 2 ^^ . Since the coincidence to accidental-coincidence ratio CAR = — is a directly 

^1^2(3) F 1 F 2(3) J 

measurable quantity, we will use CAR to quantify the temporal correlation of photon pairs 
in the following description. 

Defining the normalized detection rate of photon pairs as P pa i r = Pi£s/(ViVi)i an d using 
Eqs (1181) . ( 1191) and ( 1221) . we obtain the expression of the coincidence to accidental-coincidence 
ratio: 

CAR={ 9 ^M) T ^ + p ^% + a?) , (24) 

where a' s ^ = a s (i)/a p . It is clear that CAR decreases with the increase of P pa i r , which 
determines the brightness of HSPS. To further demonstrate the dependence of CAR, we 
then plot CAR as a function of the bandwidth of signal and idler photons. As shown in Fig. 
2, at a certain rate of P pa i r , CAR obviously increases with the increase of &' s uy Looking 
at Fig. 2 more carefully, we find that for the idler (signal) photons shaped by cr[^, if the 
bandwidth of signal (idler) photon a' s(i) is less than /2 + er /2 s, CAR increases with the 



increase of g' s u\] otherwise, CAR decreases with the increase of &' s uy Therefore, the relation 
described by Eq.( l2"4"|) is the basis for optimizing the quantum correlation between signal and 
idler photon pairs, which is the key to reliably generate heralded single photons. 

To obtain the second order intensity correlation function of the HSPS, it is necessary to 
know the theoretical expression of the triple coincidence probability between SPD1, SPD2 
and SPD3 in the same time slot: 



P 123 (0) = / dhdhdh^E^^E^^E^t^E+it^E+^Etit^^) 

dti(0| |0> / dt 2 dt 3 (0\E 2 (t 2 )E^(t 3 )E+(t 3 )E+(t 2 )\0) 



+ 2 J dt 2 (0\E 2 (t 2 )E+(t 2 )\0) J dM*3|(0|£r(*i)# 3 (*3)|0)| 2 

+ 2 J rftidtadtg <0J £7£-(*i)^(*a) l0> <0j ^(^^(ta) j0> <0| ^+(t 3 )^(*i) |0> 

= W3 + ^ 1 ^ + (sf ) -l)(W3 + ^ 1 ^) I (25) 
where g^ = 1 H — , 1 , a can be considered as the normalized second order intensity cor- 



v / i+-_ 

relation function of individual signal field 241, and the coefficient f' = — f=^== means 



10 



the collection efficiency to the two pairs events. One sees that triple coincidence probability 
P 123(0) is a sum of the three terms: the first term is the accidental coincidence between 
SPD1, SPD2 and SPD3; the second term stems from the coincidence between a signal pho- 
ton and a photon pairs; and the third term is associated with the bunching effect in the 
heralded signal band. If the bandwidth of the heralded signal a' s is very broad, the sum of 
the first two terms would be a good proximation of ^123(0); if not, however, the bunching 
effect must be taken into account. 

Having obtained the expression of the various counting probabilities of the SPDs, we 
are able to compute the second order intensity correlation for the heralded signal field, 
conditioned on observing an idler photon count in the time slot t k, i.e., 

9c (U) -p 13 (o)p 12 (o)- w 

Substituting Eqs. ( IT8l) . ( 1221) and (1251) into Eq.()26l). and using the approximation -pj + 



25j . we have 



_____ 

Eq. f )27|) shows that the value of g c (0) is determined by the second order correlation 
function of the individual heralded signal field and quantum temporal correlation of the 
photon pairs. 

Although the condition gi 2 \o) ~ is desirable for a single photon source, the desirability 
can not be satisfied unless CAR ^> 1. To achieve CAR ^> 1, according to Eq. ([21]), one 
should either decrease the brightness characterized by P pa i r , or increase a' s and a[. Because 

a decreased brightness is adverse for the applications of HSPS, to figure out how to optimize 

(2) 

the photon statistics, we plot g c (0) as a function of the bandwidth of signal and idler 
photons for P pa i r = 0.02 pairs/pulse. As shown in Fig. 3, the value of gc (0) decrease with 
the increase of a' s and and the asymmetry, owing to the bunching effect of heralded signal 
photons, can be observed in the pattern of the contour-plot. 



B. Heralding efficiency 



In addition to the intensity correlation for the heralded signal field, another important 
parameter to characterize a HSPS is the heralding efficiency, H . Given the detection of one 
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member of a photon pair, H is the probability that the twin photon is actually present in the 
output of HNLF, which can be obtained by measuring the conditional detection efficiency 
of heralded photons, and then correcting for the losses in the heralded-photon analysis zone 
(the shaded area shown in Fig. 1). If the conditional detection efficiency i]d is simply defined 
as the ratio between true coincidence and the trigger photon detection rate, i.e. 

A 2 (3)(0) - AAp) i * (0 „, 

Vd = — 5 — = r%%(3)&, (28) 

the heralding efficiency can be expressed as 



a. 



H = £ s = . (29) 
'2 + a' 2 + a' 2 



Eq. ( 129]) shows the heralding efficiency H is equivalent to the collection efficiency of the 
photon pairs for the heralding trigger photons detected in idler channel 23(. Figure 4 plots 
H as a function of the bandwidth of signal and idler photons. One sees that the broader 
the bandwidth of F2 in signal field relative to that of the filters in pump and idler fields, 
the higher the heralding efficiency. 



C. Indistinguishability 

Apart from the photon statistics and heralding efficiency, indistinguishability of HSPS 
is important as well. Because of the independent nature of the photons, in order to obtain 
high visibility in quantum interference between multiple sources, photons must be in iden- 
tical single mode so that they are indistinguishable. Although the spontaneous parametric 
emission processes are independent in different sources, the timing provided by the ultra- 
short pump pulse can be used to realize the required synchronization. For the photon pairs 
generated in fiber based sources, spatial overlap is ensured, spectral overlap can be achieved 
by the use of similar filters, and polarization can be matched by a polarization controller. 
However, the spectral of signal and idler photon are usually correlated, and the temporal 



mode structure is messed up [24|, |26( . To preserve the temporal indistinguishability, optical 
filtering is necessary so that the detected signal and idler photons are spectrally factorable. 

According to the field operators of signal and idler photons shown in Eqs. (II lj) and (112p . 
photon pairs in spectral factorable state can be obtained by using two methods. One method 
is applying a narrow band filter (NBF) in idler band, so that the heralding photons can be 
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viewed as in a single mode. For example, when bandwidth of the filter in heralding idler band 
is a[ = 0.3, the second order coherence function of the individual field gf^ — (1 H — , 1 /2 ) 



is about 1.98 24|, indicating the temporal mode of the individual idler photons can be 
considered as a single mode j^. Under this condition, no matter what kind of filter is 
applied in signal band, the HSPS is indistinguishable. The other method is inserting a NBF 
in signal band, so that the heralded photons is in a single mode. For example, when the 
filter in heralded signal band is a' s = 0.3, the HSPS is indistinguishable, with no relevant 
to the filter in idler band. Substituting the source parameters given by the two examples 

_____ /2"\ 

into Eqs. ( |24|) . ( |27|) and ( |29l) . we calculate the corresponding g c (0) and H for the pair rate 
Ppair = 0.005 pairs/pulse. As shown in Fig. ((5]), in the sense of evaluating the source by 
looking at the value of gc (0) and H, indistinguishable HSPS obtained by using NBF in the 
heralding band is slighter better, because in this case, the contribution of the bunching effect 
of the heralded signal photons (see Eq. (25)) can be made smaller, and heralding efficiency 
can be adjusted to be higher by broadening the bandwidth of heralded signal photons. 

It is well known that indistinguishable HSPS realized by using NBFs will severely limit the 
brightness of single photons. However, we would mention that the even worse characteristic 
of the HSPS is the increased value of gi 2 \o)- Comparing Fig. 5 with Fig. 3, we find 
that although the brightness in Fig. 5 is only 1/4 of that in Fig. 3, the smallest g¥\o) 
achievable in Fig. 5 is higher than that in Fig. 3. This is because the bandwidth setting 
of signal and idler photon for a indistinguishable HSPS are not optimized for obtaining a 
maximized CAR. For the individual signal and idler photons, the portion that their twin 
photons are not within the bandwidth of idler (signal) filter is not minimized. Moreover, the 
indistinguishable HSPS achieved by using a NBF in heralded band, will additionally suffer 
from the low heralding efficiency. Therefore, to develop a indistinguishable heralded single 
photons with high quality, direct generation of photon pairs in spectral factorable state is 



desirable 



12|. 



Because we have previously verified the indistinguishability of a fiber based HSPS by using 
NBF in heralded band and experimentally st udy ing the forth-order interference between 



independent sources, as demonstrated in Ref. |24j and 
the experiments described in the next section. 



we will not discuss this issue in 
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D. Interplay of photon statistics, heralding efficiency and indistinguishability 

From the above analysis of the dependence of photon statistics, heralding efficiency and 
indistinguishability, one sees the three criteria of the HSPS can not be optimized simulta- 
neously. For a HSPS with a certain brightness, i. e., the pair rate P pa i r is kept at a certain 
level, according to Eqs. (|27|) and ( 129|) . enlarging the bandwidth of signal and idler photons 

(2) 

helps to reduced gc (0) and increase H, but with a cost of decreased indistinguishability. 
On the other hand, to obtain a good indistinguishability, a NBF in signal or idler band is 
required. This will result in a relatively increased gi 2 \o), although H with a higher value 
can be maintained by optimizing the source parameters. 

It is worth noting that the brightness of HSPS, characterized by P pa i r , is also one of the 
concerns for exploring its practical applications. There are two approaches to make the 
source brighter: one is using signal and idler photons with a broader bandwidth (see Eqs, 
f TT8]) -f l22]) ). the other is to increase the gain of SFWM by increasing the pump power and 
the length of fiber. The former will result in a HSPS with a relatively low gc 2 \o) and higher 
H. However, in addition to a poor indistinguishability, this kind of source is not suitable for 
studying the Fock state with the photon number greater than 1, since the gain of SFWM 
(per mode) is low, and the two-pairs events is relatively rare. While the latter will lead to a 
increased 5^(0). But this kind of source is useful for the investigations involving the multi- 
photon pairs events. Moreover, to describe the latter more precisely, the field operators of 
signal and idler photons shown in Eq. (11) should take a higher order approximation. 

IV. EXPERIMENTS AND RESULTS 

A schematic of our experimental setup is shown in Fig. 6. Signal and idler photon pairs 
at wavelengths of 1544.53 nm and 1531.9 nm, respectively, are produced by launching a 
pulsed pump centering at 1538.9 nm along the fast or slow axis of the 20-meter-long HNLF, 
whose zero-dispersion wavelength, dispersion slope, and birefringence (An) are measured 
to be about 1525 ± 5 nm, 0.02 ps/nm 2 -km, and 2.5 x 10 -6 , respectively. The fiber po- 
larization controller (FPC1) is used to adjust the polarization of pump. The HNLF with 
7 = HW~ l km~ 1 is immerged in liquid nitrogen to suppress Raman scattering. FPC2 and 
fiber polarization beam splitter (FPBS) placed at the output port of HNLF are used to select 
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the co-polarized photon pairs and to reject the residual Raman scattering which is cross- 
polarized. The loss of HNLF mainly comes from fusing loss ( 0.8dB) due to core diameter 
mismatch between HNLF fiber and spliced standard single-mode fiber. 

The pump with a 10-ps-duration and a repetition rate of about 41 MHz is spectrally 
carved out from a mode-locked femto-second fiber laser. To achieve the required power, the 
pump pulses are amplified by an erbium-doped-fiber amplifier (EDFA), photons at the signal 
and idler wavelength from the fiber laser that leak through the spectral-dispersion optics 
and from the amplified spontaneous emission from the EDFA are suppressed by passing the 
pump through a tunable filter Fl with a FWHM of about 0.3 nm. 

The dual-band filter F2, providing a pump-rejection ratio in excess of 120 dB, is obtained 
by cascading four WDM filters (Santec, model WDM-15). To investigate the influence of the 
bandwidth of detected signal and idler photons, F2 is realized by using different combination 
of 100 GHz and 200 GHz WDM filters, so that the FWHM in signal and idler bands could be 
either 0.6 and 0.6 nm, or 1.1 and 0.6 nm, or 1.1 and 1.1 nm, respectively. For convenience, 
we label the three kinds of F2 as F2j, F2 n and F2 ni) whose spectra are plotted in inset of 
Fig. 6. 

The quantum efficiencies of SPD1, SPD2 and SPD3 are about 20%, 12% and 17%, re- 
spectively, whose corresponding dark-count probabilities are about 1.9 * 10~ 5 , 2.1 * 10~ 5 
and 5.9 * 10~ 5 counts/pulse, respectively. The electrical signals produced by the SPDs in 
response to the incoming photons (and dark counts) are acquired by a photon-counting 
system. Thus, single counts in both the signal and idler bands and coincidences acquired 
from different time slots can be recorded. To determine the total detection efficiencies for 
idler and signal channels, besides the efficiency of the corresponding SPDs, the efficiencies of 
the HNLF 70%, the filter F2 (including F2j, F2jj and F2m) in signal and idler channels, 
~ 24% and ~ 52%, 50/50 fiber coupler, and other transmission components about ~ 90% 
should be taken into account. 

Before characterizing the HSPS, we first describe the quantum correlation of photon pairs 
generated in HNLF by using SPD1 and SPD2 to record single counts and coincidence rate 
when dual band filter is set to be F2j, F2jj and F2m, respectively. Figure 7(a) shows the 
number of scattered photons in idler band per pump pulse, iVj as a function of the average 
pump power, P ave , and the measured data is fitted with the iVj = siP ave + S2-PiL, where si 
and s 2 are the linear and quadratic coefficients, which respectively determine the strengths 
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of RS and FWM in HNLF. Figure 7(b) shows the measured coincidence and accidental 
coincidence rates of signal and idler photons, produced by same pump pulse and adjacent 
pump pulses, respectively, as a function of the average pump power. It is clear, the quadratic 
part s 2 Pa Ve in Fig. 7(a) is equivalent to Pi in Eq. (18). Thus, for each kind of dual band 
filter, we can deduce the rate of signal and idler photons via SFWM, the normalized detection 
rate of photon pairs, P pa i r , and the ratio between the coincidence and accidental coincidence 
rates, CAR. Figure 7(c) shows the value of CAR versus Ppair after the effect of Raman 
scattering is subtract. One sees that for a certain value of P pa i r , CAR increases with the 
increase of the bandwidth in the detected signal and idler fields. This is consistent with the 
theory predictions in Fig. 2. 

We then use the detection signal of SPD1 to herald the existence of single photons in signal 
band, and measure its conditional second order correlation function g c (0) by recording 
the single counts and coincidences between SPD1, SPD2, and SPD3 at various conditions. 
When the detected photons of the three SPDs originate from the same pump pulse, we 
record the triple coincidence at different pump power levels by using F2j, Flu and F2m, 
respectively. During the measurement, the single counts of SPD1, SPD2, and SPD3, and 
two-fold coincidences between two arbitrary SPDs are also recorded. The corresponding 
value of g c (0) can be obtained according to Eq. (1261) . To figure out the dependence of 

(2) 

g c (0) upon the brightness in different cases, we extract the normalized detection rate of 
photon pairs, P pa i r from the measured single counts and two-fold coincidences, and sketch 

(2) (2) 

g c (0) as a function of P pa i r - As shown in Fig. 8 (a), for a certain P pa i r , 9c (0) obtained 
by using F2jjj, with FHWM of 1.1 nm in both signal and idler channels tends to be the 

(2) 

smallest, while g c (0) obtained by using F2j, with FHWM of 0.6 nm in both signal and 
idler channels, tends to be the biggest. The result is qualitatively in accordance with the 
theory expectation in Fig. 3. We note that with the reduction of P pa i r , the discrepancy 
between the three sets of data become unapparent. This is because the contamination of 
Raman scattering is getting more severe with the reduction of the pump power and with 
the broadening of the bandwidth of signal and idler photons, as illustrated by the fitting 
parameters and curves in Fig. 7(a). Therefore, correcting the measured results of g¥\o) by 

(2) 

subtracting the influence of Raman scattering, we plot the g c (0) as a function of P pa i r , and 
compare the modified data with the calculated results. The calculation is carried out by 
substituting the source parameters into Eqs. (12^1) and (1271) . in which no fitting parameter 
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is used, and the FWHMs of the Gaussian shaped filters are chosen to be the same as that 

(2) 

of the Super Gaussian shaped F2. As shown in Fig. 8 (b), the corrected g^ (0) agrees with 
calculated results, although there is departure between them, which is originate from the 
differentia of the filters used in experiments and theoretical model, respectively. 

Finally, we characterize the heralding efficiency of the HNLF-based HSPS. Using the data 
in Fig. 7(b), we deduce the true coincidence, Rt rue , which is the difference of coincidence 
and accidental coincidence rate between SPD1 and SPD2. To obtain the heralding efficiency 
H, according to Eqs. (|28|) - (|29|) . we first divide Rtrue by the single counts of SPD1, which is 
the trigger photon detection rate, and then correct the result by using the total detection 
efficiency in the signal channel connected with SPD2. Figure 9(a) plots the measured if as a 
function of pump power for signal and idler photons with different bandwidth combinations. 
One sees at a certain power, the biggest and smallest value of H can be achieved by using 
signal and idler photons shaped by F2jj and F2j, respectively. In the sense of comparing 
the relative value of H obtained in different cases, the results qualitatively agrees with Eq. 
(|29l) . However, instead of being irrelevant to pump power, it is obvious that for signal 
and idler photons with a certain spectra, H increases with the increase of pump power. 
The departure is caused by the contamination of Raman scattering. To confirm this point, 
using the fitting results in Fig. 7(a), we modified the results by subtracting the influence 
of Raman scattering and compare them the calculated results obtained by substituting the 
source parameters into Eq. (!29l . As shown in Fig. 9(b), for the signal and idler photons 
with a specified spectra, the modified H does not depend on the pump power, and agree 
with theory expectations. 



V. CONCLUSION 



In conclusion, using the photon pairs produced in 20-meter long HNLF by a pulsed pump 
in 1550 nm band, we have studied the HSPS from the theoretical and experimental aspects. 
Based on the Bogoliubov transformation for a broadband pulsed pump, we derive the an- 
alytic expressions of the second-order coherence function (0) and heralding efficiency H 

(2) 

of the HSPS for the first time. The calculation reveals the dependence of g c (0), H and in- 
distinguishability upon the source parameter, and the validity of the theoretical calculations 
are proved by the experimental results. We believe the theoretical frame used to model the 
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fiber based HSPS is also useful for describing other spontaneous parametric process pumped 
with a mode-locked laser. 

The HSPS presented in this paper is an all-fiber source with compact size and freedom 
from misalignment. Although heralded single photons in 1550 nm telecom band are desir- 
able for quantum communication, we would like to mention that limited by the detection 
techniques, heralding photons in visible band are more practical at the current stage 11]. 
However, the use of 1550 nm heralding photons brings us a lot of conveniences in inves- 

(2) 

tigating the dependence of g c (0) and H under different experimental conditions, because 
various kinds of off-the-shelf fiber components with high quality and low cost are commer- 
cially available in this band. 

We think the HSPS can be improved from the following aspects. Firstly, the discrepan- 
cies between the directly measured experimental results and theory calculations due to the 
existence of Raman scattering can be reduced or wiped off by cooling the HNLF to 4 K 29] 
or by using photonic crystal fiber based photon pairs with large detuning to further mitigate 



Raman effect 



9-11 



22 ] . Secondly, the current transmission loss experienced by heralded 



signal photons can be reduced to less than 1 dB by minimizing splicing loss and using high 
quality WDM or fiber-bragg-grating filters. Thirdly, the emission rate of the heralded single 
photons can be be dramatically increased by pumping the HNLF with a mode-locked fiber 
laser operating at a repetition rate above 10 GHz and by using SPDs with high speed and 
high efficiency 30, [3l| . Finally, to further increase the single photon rate while maintaining 

(2) 

a low value g c (0), similar to multiplexed SPDC scheme used to better approximate a single 
photon source on demand [32j, it is possible to operate an array of simultaneously pumped 
SFWM sources. In this case, in the sense of obtaining a reduced g c (0) and increased H 
at a certain pair rate P pa ir for each SFWM source, our study on characterizing the HNLF 
based HSPS would be useful for optimizing parameters. 
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FIG. 1: (Color online) Conceptual representation of the scheme for generating heralded single 
photons. F2, dual-band filter. 




FIG. 2: (Color online) Contour-plot of quantum correlation of photon pairs CAR as a function 
of the bandwidth of signal and idler photons. The calculation is carried out at the pair rate 
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FIG. 3: (Color online) Contour-plot of g ( c '{0) as a function of the bandwidth of signal and idler 
photons for P pa ir =0.02 pairs/pulse. 
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FIG. 4: (Color online) Contour-plot of heralding efficiency of the heralded single photons in signal 
field H as a function of the bandwidth of signal and idler photons. 
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FIG. 5: (Color online) §c (0) and H of indistinguishable heralded single photons with a pair rate of 
Ppair = 0.005 pairs/pulse for (a) heralding photons is in a single mode, realized by setting <r- = 0.3, 
and (b) heralded photons is in a single mode, achieved by setting a' s = 0.3. 




FIG. 6: (Color online) Schematic of the experimental setup. The inset is the spectra of dual band 
filter F2, stars and diamonds represent the spectra with FWHM of 0.6 and 1.1 nm, respectively. 
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FIG. 7: (Color online) (a) The number of scattered photons in idler band per pulse, iVj, as a function 
of the average pump power, P ave - Thick and thin solid curves are the fitting of iVj = siP ave +S2Pa Ve , 
with si = 0.061, s 2 = 0.027, and si = 0.030, s 2 = 0.012, for F2 in idler band with FWHM of 1.1 and 
0.6 nm, respectively, (b) Coincidences rates (represented by solid circles, triangles and diamonds) 
and accidental coincidence rate (represented by hollow circles, triangles and diamonds) of the 
detected signal and idler photons versus P aV e- (c) The deduced ratio between the coincidence and 
accidental coincidence rates, CAR, as a function of the normalized photon pair detection rate, 
Ppair, after the influence of Raman scattering is excluded. The solid curves are only for guiding 
eyes. 
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FIG. 8: (Color online) (a) Measured and (b) corrected g c (0) of the HNLF-based HSPS versus 
brightness, characterized by the normalized photon pair detection rate P pa ir- The size of the error 
bar of Ppair is determined by the uncertainty of the total detection efficiencies, and the uncertainty 

(2) 

of (0) lies on the statistical fluctuation of photon counting. The theory curves in (b) are 
calculated by using Gassian shaped filters with FWHM the same as that of F2. 
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FIG. 9: (Color online) (a)Heralding efficiency of HNLF-based HSPS, H, as a function of average 
pump power, (b) The corrected heralding efficiency versus the average pump power when the 
Raman photons in idler bands are subtracted. The straight lines in (b) are the theoretically 
calculated results, in which the FWHMs of Gaussian shaped filters are the same as that of F2. 
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